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Abstract. Recursion relations are used to exactly calculate the partition function of 
a canonical ensemble in which all additive charges as well as the total isospin are 
strictly conserved. The ensemble can consist of particles that obey either classical 
or quantum statistics. Recursion relations are also employed to compute observables 
such as multiplicity distributions and isospin fluctuations. 
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1. Introduction 

Results from a canonical ensemble that obeys exact conservation of additive charges and 
isospin can be useful in a variety of applications. For example, the unexpected isospin im- 
balances of CENTAURO events [ [l| could be explained by employing a Gaussian thermal 
source and restricting the pion emission to isoscalar pairs only. The pion multiplicity distri- 
butions become broader as compared to those of independent pion emission, thus resulting 
in larger fluctuations [ ph. This finding could be further substantiated by replacing the 
Gaussian thermal source by a canonical ensemble and the ad-hoc assumption of isoscalar 
pairs emission by conservation of total isospin and its projection. 

Another potential area of application is the modeling of relativistic heavy ion colli- 
sions. Here, hydrodynamical models are used to describe the early stages of the collision, 
when the density is so high that the particles behave collectively as a fluid. As the collision 
progresses the density falls, the system breaks up, and hydrodynamics is not applicable 
anymore. Therefore, the break-up and later stages of the collision have to be described by 
other means such as, for example, microscopic transport models [ EJ. When switching to 
a microscopic transport model from a macroscopic hydrodynamical prescription, particles 
must be formed from the fluid such that momentum and energy as well as isospin and all 



ISBN 963 00 0000 
© 2002 EP Systema, Debrecen 



2 



Sen Cheng 



other charges are conserved. Currently, most schemes for generating particles are based on 
the grand canonical ensemble. However, exact conservation within each event and within 
each computational cell is essential for studies on fluctuations and balance functions [ph. 

Partition functions for canonical ensembles conserving charges can be written quite 
easily, but they are very difficult to compute for large systems because the number of par- 
titions that have to be summed over increases rapidly [ n]. This difficulty can be overcome 
by rewriting the partition function as a recursion relation. These recursion relations were 
applied to ensembles of classical particles in [ and to Fermi systems in [ |[|. 

In this work we add total isospin conservation to the canonical ensemble and derive 
the appropriate recursion relations. Furthermore, isospin fluctuations in the canonical en- 
sembles, which, for example, are relevant for CENTAURO event studies, are calculated. A 
Monte Carlo algorithm that employs these recursion relations to generate particles is part 
of ongoing work and will be presented in the future. As our focus is the development of 
the computational technique, we show sample calculations for simple systems rather than 
comparisons to experimental results. 



2. Recursion Relations with Classical Statistics 



2.1. Partition Functions 

The partition function for a canonical ensemble of A particles conserving a charge Q, where 
Q could be a vector containing several charges, is defined as 

Z A . e = £(a| e -P», (1) 

a 

where a labels some quantum state with the appropriate quantum numbers, p is the inverse 
temperature, and H is the Hamiltonian of the system. Using the single particle partition 
function 

co* = £exp(-e[V), (2) 

i 

which is simply a sum of Boltzmann weights over all available states i, the partition func- 
tion can be written in a more tractable form as a sum over partitions subject to particle 
number and charge conservation 

N 

z ^ = I n£r (3) 

N is the number of particle types, V* is the number of particles of type k, and is the charge 
of one particle. The particle number indicates how many times a particle contributes to 
the main conserved quantity, e.g., in the case of multi-fragmentation, would be the mass 
number of the nuclei. 

The computational challenge of Eq. (Q) is hidden in the number of partitions to be 
summed over. Assuming = 1 for all k, the number of partitions is 
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which is prohibitive to compute for large A or N. However, summing over all partitions can 
be bypassed when using recursion relations as derived in [ 

Za -Q = J Y, a i (S) } ZA - a i>Q-ir ( 5 ) 

If intermediate values are stored, Eq. (|j) requires a linear runtime in both, A and N. 

To conserve isospin as well, a sum over all possible isospin configurations for a given 
partition, {v^}, and isospin weights W (I,M\ {Vk}) have to be included in the partition func- 
tion 

N ^ 

£W,e = L L^(/,M|{va)n^t- (6) 

ljy k a k =A-Xy k q k =Q){v k } k=\ Wk ~ 

This partition function can be converted to a recursion relation by inserting jY!j=i a j v j = 1 



1 N 0)^ 1 fo v * 

n A j,M,Q = T I aj*j L E W (/,M|{v*» n -fr 

J=l <£v*a*=^;£v Mjt =Q){v t } ^ 1 ^t^y v *- 

The summation indices can be switched 



(7) 



y' k = { Vk , (8) 

* \ Vfc-l = ; 
if the sum over isospin configurations is split: 

ij+l 

X>(/,M|{v t }) = £ (Ijm j ;l' ) M-mj\IM) 2 Y J W(l',M\{v' k }). (9) 

K} /'=|M- m; | {v^} 

Finally, the recursion relation with conserved total isospin turns out to be 



1 I+Ij 

&A,I,M,Q = tL^' 03 ^' E ( I j m j'> I '> M - m j\ IM ) 2 £ 1 A-a j ,I>,M-mj,Q-q j - 

j I' = \M-mj\ 



(10) 



2.2. Multiplicity Distributions 

We derive a recursion relation for the multiplicity distribution of a set B of particle types. 
The probability of observing n particles of a type contained in the set B can be expressed 
in terms of factorial moments [ [7|] 

i>„ \ l n).n. 

for which recursion relations can be found more easily. Following the procedure in the 
previous subsection the recursion relation for the factorial moment is found to be 

F B A,i,M(n) = £ (d k Y,(hm k ;l',M -m k \IM) l Q ' A ~°*> I '> M ~ mk F BtA _ ak/jM -m k (n - 1), (12) 

"A. ' ' ' 
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which leads to the recursion relation for the multiplicity distribution 

PB.A.i.M{n) = - £ (O k Y(I k m k ;l',M-m k \IM) 2 A ~^^- m ^ P BA _ a f M _ m (n-\). (13) 

To test Eq. ( |l3| ) we implemented a program that included the pseudo-scalar and vector 
mesons as well as their decays. Figure [l] shows the multiplicity distribution for neutral 
pions: the system was constrained to be in an isospin singlet at a temperature of 170 MeV 
and had to produce a total of 20 pions. For this and all following sample calculations we 
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Fig. 1. Multiplicity distribution of neutral pions for several volumes. The calculation 
includes all pseudo-scalar and vector mesons and their weak decays. The system was 
constrained to an isospin singlet at a temperature of 170 MeV. The total number of pions 
was a fixed A = 20. 

need to make an assumption about the energy states that are available to the particles and 
that have to be summed over in the single particle partition function of Eq. (Q). We required 
that the linear momentum k in one dimension satisfy the condition 

kR = %(n + ^j, (14) 

where n = 0, 1,2, . . . and R is the size of the system. The half integers on the RHS, instead 
of the more usual integers, were inspired by the MIT bag model and lead to a reduction of 
zero-point surface energy effects. This becomes important for systems that are confined to 
a small volume. By assuming that our model systems are confined to a cube of volume V, 
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we obtain the energy states 



£n,m,l = y |j [(" + V2) 2 + (m + 1/2)2 + (, + 1/2) 2] +M 2 7 (15) 

where n,w,Z = 0, 1,2, . . ., M is the mass of the particle, and R = V 1 / 3 . 

The even-odd asymmetry in Fig. [j]for small volumes that vanishes for larger volumes 
can be understood in the following way. At small volumes and low temperatures the system 
minimizes the number, not the mass, of particles produced; it prefers to produce heavy 
particles that decay into multiple pions rather than producing pions or kaons. These heavy 
particles (co, T), T|') decay into an odd number of pions (and also into an odd number of 
neutral pions). Because the number of heavy particles has to be even to make 20 pions, as 
we imposed, the number of neutral pions is strongly favored to be even. At higher volumes 
or temperatures low mass particles dominate and the even-odd asymmetry vanishes. 
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Fig. 2. Total number density of a system of only pions in an isosinglet at a temperature 
of 170 MeV. Shown are calculations for a grand-canonical and a canonical ensemble with 
classical statistics. 

The factorial moment from Eq. ( |l2| ) corresponds to the density, when n = 1. For 
a system of pions in an isosinglet, particle production is significantly suppressed at low 
volume, but converges towards the grand-canonical result at large volumes (Fig. 0). 
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3. Recursion Relations with Quantum Statistics 

For a canonical ensemble of quantum particles with charge conservation the recursion re- 
lation for the partition function was given in [ ||] : 

2 A.Q = H ^A-n,Q-q v (16) 

A k n=\ 

where C« = (±1)" _1 £;exp(— nE, /T), To derive a partition function Qa.i.m.Q ma t con- 
serves total isospin as well, we first observe that, if the Hamiltonian of the system is inde- 
pendent of M, then so is the partition function. Then the partition function conserving total 
isospin can be related to one that conserves only the projection Z^.m.q = Hi>m^a.i.m,q- 
Therefore 

&A,I,M,Q — Zam.q — Zam+i.q- (17) 
To calculate densities the 2-point function is needed 

(ajaj) = ^^(-l)"- 1 exp(-«e,/r)Z A _„, (18) 

A n 

which, however, can only be used when conserving additive quantities, but not isospin. The 
4-point function serves to calculate higher moments 

{a]a]a k a,) = ^-(S a 6 jl ±S tt 5 J/ ) £ (-l)"^exp(-n i e / /T-n^/r)Z A _„ i _„ r (19) 

L A njMj 

4. Isospin Fluctuations 

For classical particles the isospin fluctuation 

a=£</|(» + + «_-2«o) 2 |/), (20) 

7=0 

where n+, n_, no are the densities of the respective pions, can be readily computed with the 
recursion relations derived above. However, for quantum particles we cannot compute the 
density, let alone the higher moments, when conserving total isospin. Therefore, we need 
to rewrite Eq. ( p0| ) in terms of isospin projection states. 

The operator in Eq. (p0|), itself a product of spherical tensor components, can be 
written as a decomposition into other spherical tensor components A jm 

(«++«_ -2« ) 2 = 672o72o = 6^(20;20|7M)A iM . (21) 

J,M 

By the Wigner-Eckart theorem only Aoo contributes when contracted between isosinglet 
states, a = 6£/ = o(/|Aoo|/), which, after some algebra, becomes 
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Fig. 3. Isospin fluctuations as a function of volume for a system of pions at a temperature 
of 170 MeV with a fixed number of pions, A = 20. Shown are calculations with classi- 
cal and Bose-Einstein statistics. The isospin fluctuations are divided by 2A, the expected 
isospin fluctuation for a gas of independent pions. 

In Figs. H and ^ we compare classical to Bose-Einstein statistics in calculations of 
isospin fluctuations for a system of pions at T = 170 MeV. The fluctuations are divided by 
2A, the expected fluctuation for an independent pion gas, to highlight the difference. When 
the system is fixed to contain 20 pions (Fig. [|), the isospin fluctuation for classical particles 
does not depend on the volume, as the relative composition of the classical pion gas does 
not. Bose-Einstein particles, however, prefer to be all in the same state, when the volume 
becomes small, thus enhancing the isospin fluctuations. When the number of pions is not 
fixed, the isospin fluctuations converge towards the value for an independent pion gas as 
the system increases in size (Fig. Q). 

5. Conclusions 

Recursion relations have been derived for a canonical ensemble in which all additive charges 
as well as the total isospin are strictly conserved. In the case of classical statistics we de- 
rived recursion relations for the partition function and multiplicity distribution. For quan- 
tum statistics we obtained a recursion relation for the partition function and density. Work 
is under way to derive the recursion relation for the multiplicity distribution in the latter 
case. 

These recursion relations can be applied to an arbitrary number of particle types and 
conserved additive charges. They can be computed correctly within the numerical precision 
of the computer in a reasonable runtime as our sample calculations show. We can now 
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Fig. 4. Isospin fluctuations as a function average pion number for a system of pions at a 
temperature of 170 MeV. Shown are calculations with classical and Bose-Einstein statistics. 
The isospin fluctuations are divided by 2(A), the expected isospin fluctuation for a gas of 
independent pions. 

move forward and apply this technique to physical problems and in Monte Carlo algorithms 
that generate particles with exact conservation of additive charges and isospin. We also 
derived a method to quickly calculate isospin fluctuations, which makes use of the recursion 
relations we obtained. 
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